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ABSTRACT 


This  report  provides  a  summary  of  an  extended  program  of 
theoretical  and  experimental  research  on  turbulent  boundary 
layers  under  supersonic  free  stream  conditions.  The  program 
began  with  Investigations  of  boundary  layer  behavior  at 
various  Mach  Numbers  for  smooth,  adiabatic  surfaces.  Later 
phases  Involved  the  addition  of  various  types  of  surface 
roughness,  with  consideration  in  the  final  aspects  being 
given  to  the  combined  effects  of  roughness  and  heat  transfer. 
The  experimental  program  involved  the  utilization  of  velocity 
profile  measurements  and  the  momentum  deficit  method  of  calcu¬ 
lating  skin  friction,  the  application  of  both  the  floating 
element  skin  friction  balance  and  the  Preston  tube  for 
deteminations  of  local  shear  stress  values,  and  the  utiliza¬ 
tion  of  plug-type  calorimeters  for  deteminations  of  heat 
transfer  rates. 
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PREFACE 


The  ARL  program  of  boundary  layer  research  was  Initiated  In 
1947  in  order  to  provide  basic  technical  data  on  the  flow  charac¬ 
teristics  of  turbulent  boundary  layers  under  supersonic  free  stream 
conditions.  The  Immediate  application  of  such  information  was 
contemplated  as  related  to  aerodynamic  design  problems  of  guided 
missiles  and  high  speed  aircraft,  both  of  which  were  beginning  to 
attain  operational  values  of  flight  Mach  Numbers  greater  than 
unity.  Bnphasis  was  placed  on  the  turbulent  boundary  layer  because 
it  was  felt  that  surface  conditions  on  a  typical  flight  vehicle 
would  be  most  likely  to  result  in  such  situations,  rather  than  in 
the  laminar  type  of  flow  usually  related  to  lower  ranges  of  the 
flight  Reynolds  Number  and  to  extremely  small  values  of  surface 
roughness . 

The  ARL  aerodynamic  research  program  as  a  whole  developed  out 
of  the  writer's  participating  during  19^5  in  the  Bumblebee  program 
of  the  Applied  Physics  Laboratory  of  The  Johns  Hopkins  University. 
After  completion  of  that  assignment  and  return  to  The  University 
of  Texas  at  Austin,  arrangements  were  made  for  a  continuation  of 
certain  phases  of  the  Bumblebee  program  at  the  latter  institution. 
Such  activities  were  initially  conducted  under  the  aegis  of  an 
interdisciplinary  activity  originally  designated  as  the  Defense 
Research  Laboratoiy  and  later  renamed  Applied  Research  Laboratories, 
its  present  title. 

Support  for  these  DRL/ARL  research  activities  was  provided 
initially  by  the  Navy  Bureau  of  Ordnance  through  a  series  of  research 
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contracts  which  were  monitored  by  the  Applied  Physics  Laboratory. 

A  later  reorganization  of  the  Navy  Department  resulted  in  a  transfer 
of  responsibility  for  these  programs  to  the  Naval  Air  Systems  Command, 
with  technical  responsibility  continuing  In  the  hands  of  AFL/JHU. 
Significant  contributions  were  also  made  by  the  Air  Force  Office  of 
Aerospace  Research  and  other  Air  Force  technical  divisions,  particu¬ 
larly  In  the  development  of  a  significant  supersonic  wind  tunnel 
facility  that  made  it  possible  to  conduct  experimental  Investigations 
directly  at  the  University. 

In  addition  to  providing  useful  Information  for  the  aerodynamic 
design  problems  of  high  speed  flight  vehicles,  the  research  programs 
enabled  engineering  faculty  members  to  participate  actively  In  work 
of  current  and  future  technical  Importance.  At  the  same  time  a 
substantial  number  of  graduate  students  In  engineering  were  employed 
as  research  engineers  and  assistants,  with  arrangements  being  made 
which  enabled  them  to  use  the  results  of  their  studies  as  the  basis 
for  master's  theses  and  doctoral  dissertations.  Many  of  these 
individuals  have  now  moved  on  to  higher  level  positions  In  education. 
Industrial  organizations,  or  research  and  development  laboratories, 
where  they  continue  to  make  substantial  contributions  to  technology. 

In  view  of  the  large  number  of  individuals  involved  in  all 
aspects  of  this  program,  it  Is  Impossible  to  acknowledge  their 
contributions  on  an  Individual  basis.  They  can  be  recognized  only 
as  members  of  the  various  participating  groups  Involved  in  the 
support  of  the  program,  including  the  technical  divisions  of  the 
Navy  Department,  The  Johns  Hopkins  Applied  Physics  Laboratory,  and 
the  Air  Force  Office  of  Aerospace  Research.  The  opportunity  to 
participate  in  the  exchange  of  Information  resulting  from  membership 
in  the  technical  panels  of  the  Bureau  of  Weapons  Advisory  Committee 
on  Aeroballlstics  and  the  Navy  Aeroballlstlcs  Advisory  Committee 
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was  of  Inestimable  value.  Recognition  of  the  contributions  of 
individual  research  workers  at  The  University  of  Texas  at  Austin  is 
given  by  listing  of  the  technical  reports  for  which  they  were 
responsible. 
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NOMENCLATURE 


a  -  constant  In  roughness  function  for  flat  plates 

b  -  constant  in  roughness  function  for  flat  plates 

A  -  constant  in  "law  of  the  wall"  equation 

A  -  coefficient  in  "law  of  the  wall"  equation  for  rough  plates 

r 

B  -  constant  in  "law  of  the  wall"  equation 

Br  -  coefficient  in  "law  of  the  wall"  equation  for  rough  plates 

B'  -  coefficient  in  "law  of  the  wall"  equation  for  rough  plates 

r 

c  -  constant  in  Harkness's  cubic  temperature  function 

cf  -  local  skin  friction  coefficient 

C  -  constant  in  velocity  defect  law 

C'  -  constant  in  "law  of  the  wake"  equation 

Cp  -  mean  skin  friction  coefficient  for  flat  plates 

C_  -  value  of  C_  based  on  wall  conditions 

W  * 

Cp  -  specific  heat  at  constant  pressure 

d  -  diameter  of  Preston  tube 

d  -  limiting  value  of  Preston  tube  diameter 

cr 

Df  -  frictional  force  on  flat  plate 

e  -  base  of  natural  or  Napierian  logarithms 

E  -  Eckert  Number  =  2(1-^  )/(T  -T^ 

f(tl)  -  Coles'  function  in  "law  of  the  wall"  equation 

f^(rj)  -  "law  of  the  wall"  value  of  f(ri) 

fg^)  -  "law  of  the  wake"  value  of  f ( tj) 

f(rir)  -  roughness  function 

F(y/6)  -  velocity  defect  law  function 

h  -  heat  transfer  coefficient 

k  -  proportionality  coefficient  in  mixing  length  relation 

kadm  -  admissible  sand  grain  roughness  diameter 

kr  -  mean  roughness  diameter  or  height 

k  -  average  sand  grain  roughness  diameter 

s 

K  -  constant  in  transition  equation  for  Cp 


K' 

I 

i 
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O 


S 

t 

t 

t' 

T 


u 


u. 


-  constant  in  "lav  of  the  wake"  equation 

•  mixing  length 

-  plate  length,  measured  from  leading  edge 

-  Mach  Number  function  »  ( 7-l)M^/2 

-  Mach  Number  at  any  point  in  flow 

-  Mach  Number  at  outer  edge  of  boundary  layer 

-  Nusselt  Number  =  hi/X 

-  Reynolds  Analogy  factor 

-  Prandtl  Number  =  pC ^/X 

-  Prandtl  Nvmber  at  vail 

-  Prandtl  Number  at  insulated  vail 

-  heat  transfer  rate 

-  pipe  radius 

-  temperature  recovery  factor 

-  Reynolds  Nunber  =  pu t/\i 

-  Reynolds  Number  at  edge  of  laminar  sublayer 

•  Reynolds  Number  based  on  vail  conditions 

-  value  of  friction-distance  parameter  at  edge  of  laminar 
sublayer 

-  value  of  s  at  edge  of  laminar  sublayer  for  Insulated 
plate 

-  Stanton  Number  ■  h/p 

-  vail  thickness  for  Preston  tube 

-  temperature  function 

-  temperature  function 

-  absolute  temperature 

-  temperature  at  outer  edge  of  boundary  layer 

-  stagnation  temperature 

-  vail  temperature 

-  adiabatic  vail  temperature 

-  velocity  parallel  to  vail  surface  at  any  point  within 
boundary  layer 

-  value  of  u  at  outer  edge  of  boundary  layer 
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u  -  wall  velocity,  value  of  u  at  outer  edge  of  laminar 

sublayer 

u*  -  friction  velocity 

w(y/  )  -  Coles’  wake  function 

x  -  distance  along  plate  surface  measured  from  leading  edge 

X  -  function  of  Mach  Number  and  skin  friction  coefficient 

X’  -  function  of  Mach  Number 

y  -  distance  across  boundary  layer  measured  from  wall 

Y  -  function  of  Mach  Number  and  skin  friction  coefficient 

Y*  -  function  of  Mach  Number 


xiii 


Greek  Symbols 


0  -  exponent  In  "law  of  the  wall"  equation  for  compressible 

flow 

7  -  ratio  of  specific  heats  at  constant  pressure  and  constant 

volune 

5  -  boundary  layer  thickness 

-  momentum  thickness 

5  -  laminar  sublayer  thickness 

X  -  function  of  Mach  Number  and  temperature  ratio 

X  -  coefficient  of  thermal  conductivity 

u  -  absolute  viscosity  of  fluid 

-  value  of  p  at  outer  edge  of  boundary  layer 

v  -  kinematic  viscosity  of  fluid 

-  value  of  v  at  outer  edge  of  boundary  layer 

v  -  value  of  v  at  wall 

w 

n  -  wake  factor  in  Coles’  "law  of  the  wake" 

p  -  mass  density  of  fluid 

o  -  function  of  Mach  Number 

t  -  shear  stress 

t  -  shear  stress  at  wall 

o 

cp  -  ratio  of  flow  velocity  to  friction  velocity  =  u/u* 

qp^  -  value  of  cp  at  outer  edge  of  boundary  layer 

♦  -  function  of  a  and  X 

cd  -  exponent  in  viscosity-temperature  relation 

Subscripts 

1  -  outer  edge  of  boundary  layer 

o  -  wall  conditions 

s  -  outer  edge  of  laminar  sublayer 

w  -  wall  conditions 
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I.  INTRODUCTION 

This  report  presents  a  summary  of  research  on  turbulent  boundary 
layers  carried  on  over  an  extended  period  at  the  Applied  Research 
Laboratories  of  The  University  of  Texas  at  Austin.  The  work  was 
supported  initially  by  the  Navy  Bureau  of  Ordnance,  and  later  by  the 
Bureau  of  Weapons  and  the  Air  Systems  Command,  through  a  series  of 
contracts  with  the  Navy  Department  and  more  recently  through 
subcontracts  with  the  Applied  Physics  Laboratory  of  The  Johns  Hopkins 
University.  Current  activities  are  covered  by  AFL/JHU  Subcontract 

In  its  early  stages  the  research  program  waB  concerned  with 
the  determination  of  skin  friction  and  boundary  layer  characteristics 
on  a  smooth  flat  plate.  The  experimental  conditions  selected  for 
study  were  supersonic  free  stream  velocities  and  adiabatic  conditions 
insofar  as  heat  transfer  was  concerned,  the  primary  objective  being 
to  provide  quantitative  data  on  skin  friction  suitable  for  applica¬ 
tion  to  missile  design  problems  over  a  moderate  range  of  Mach 
Numbers.  Turbulent  flow  conditions  were  chosen  because  it  was  felt, 
at  that  time  at  least,  that  the  boundary  layer  flow  on  any  operational 
vehicle  would  probably  be  of  this  character.  Heat  transfer  effects 
were  omitted  from  these  early  studies  in  order  to  provide  a  relatively 
simple  approach  to  these  problems.  The  initial  phase  of  the  test 
program  was  conducted  in  the  supersonic  wind  tunnel  of  the  Ordnance 
Aerophysics  Laboratory  at  Daingerfield,  a  facility  which  made  it 
possible  to  achieve  a  fairly  high  range  of  Reynolds  Numbers  as  well 
as  a  modest  range  of  low  supersonic  Mach  Numbers.  In  1957  the  QAL 
wind  tunnel  was  decommissioned  and  all  testing  was  transferred  to  a 
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new  facility  developed  at  ARL  with  support  from  the  Air  Force  Office 
of  Aerospace  Research.  The  ARL  tunnel  has  a  substantially  smaller 
test  section,  6  in.  by  7  In**  as  compared  to  the  19  in.  by  22  in. 

OAL  tunnel,  but  the  ARL  tunnel  does  cover  a  wider  range  of 
Mach  Numbers,  i.e.,  from  2.00  to  5.00,  approximately. 

During  the  period  of  operation  of  the  OAL  tunnel,  the  research 
program  was  extended  by  the  introduction  of  roughness  effects, 
initially  of  a  uniformly  distributed  character,  and  later  in  the 
form  of  transverse  Veegrooves.  Considerable  attention  was  also 
devoted  to  the  development  of  various  types  of  instrumentation  for 
the  measurement  of  local  shear  stress.  The  initial  smooth  plate 
work  had  been  accomplished  by  the  use  of  a  Pitot  probe  which  traversed 
the  boundary  layer  and  provided  a  means  of  determining  the  mean  skin 
friction  up  tcj  the  measuring  station  by  a  momentum  deficit  calculation. 
Early  analytical  studies  of  turbulent  boundary  layers  brought  forth 
the  significance  of  local  shear  stress  and  its  relation  to  the 
velocity  distribution.  This  concept  led  to  the  development  of  a 
floating  element  type  of  balance,  capable  of  measuring  the  shear 
force  on  a  small  disk  set  flush  with  the  main  boundary.  The  Pitot 
probe  method  was  also  extended  to  include  Preston’s  technique  of 
local  shear  determination  by  means  of  a  tube  resting  directly  on  the 
plate  surface.  Correlations  between  these  different  methods  were  the 
subject  of  extended  study. 

In  the  development  of  the  ARL  wind  tunnel,  one  of  the  initial 
design  considerations  involved  the  inclusion  of  an  air  preheating 
unit,  primarily  as  a  means  of  preventing  condensation  in  the  flow 
passing  through  the  test  section.  The  requirement  for  such  a 
heater  then  led  to  an  extension  of  the  capacity  of  the  unit,  so  as 
to  make  possible  the  execution  of  boundary  layer  tests  with  heat 
transfer  for  a  modest  range  of  heat  transfer  rates.  In  seme  cases 
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these  rates  were  extended  by  internal  cooling  of  the  models  being 
tested.  The  heat  transfer  tests  were  conducted  both  on  a  smooth 
plate  and  on  roughened  plates. 
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II.  THEORETICAL  BACKGROUND  FOR  INCOMPRESSIBLE  PLOW 


In  order  to  fully  understand  the  details  of  the  experimental 
programs  and  the  nature  of  the  results  obtained,  it  is  necessary  to 
provide  a  brief  outline  of  basic  boundary  layer  theory.  The  AHL 
program,  with  a  few  exceptions  which  will  be  mentioned  later,  did  not 
Include  extensive  analytical  treatments  of  boundary  layer  flows,  but 
rather  utilized  existing  methods  for  this  purpose. 

Although  many  refinements  have  been  introduced  in  boundary 
layer  analysis,  the  most  successful  theoretical  work  still  remains 
the  mixing  length  theory,  initiated  by  Prandtl,  and  later  extended 
by  Von  Kannan.  The  basic  physical  assumptions  concerning  the  flow 
are  well  described  in  any  of  the  numerous  treatises  on  fluid 
mechanics,  with  the  work  of  Schlichting  [Ref.  (l)]  representing  one 
of  the  most  comprehensive.  The  starting  point  of  the  analysis  is 
the  differential  equation  for  the  shear  stress  at  any  point  in  the 
boundary  layer,  which  Prandtl  derived  in  the  form 

t  =  -pf2(du/dy)2  .  (l) 

This  relation  is  based  on  the  concept  that  momentum  is  transferred 
from  one  layer  of  fluid  to  another,  there  being  a  gradual  retardation 
in  the  flow  velocity  as  the  boundary  surface  is  approached.  This 
basic  equation  has  been  used  in  studies  of  pipe  flow  and  for  boundary 
layer  flows  on  flat  plates  and  curved  walls,  the  latter  case  involving 
a  longitudinal  pressure  gradient. 


Preceding  page  blank 
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In  order  to  carry  out  the  integration  of  Eq.  (l),  It  is 
necessary  to  have  available  two  additional  pieces  of  infonnation. 

One  of  these  is  the  variation  in  shear  stress  across  the  boundary 
layer,  that  is,  t  as  a  function  of  y,  while  the  other  is  the  value 
of  the  mixing  length  term  l.  Prandtl  reduced  the  problem  to  a 
relatively  simple  form  by  assuming  that  in  the  immediate  vicinity 
of  the  boundary  the  shear  stress  is  constant  and  equal  to  the  value 
at  the  wall.  His  value  of  the  mixing  length  factor  was  based  on 
experiments  made  by  Nikuradse  on  flow  in  smooth  pipes  which  showed 
a  pseudoparabolic  variation  in  i  with  y,  reaching  a  maximum  at  the 
center  of  the  pipe.  Again,  focusing  his  attention  on  the  wall 
region,  Prandtl  used  an  approximate  relation  of  the  form 

I  a  ky  ,  (2) 

where  k  *  0.4  is  the  slope  of  the  curve  of  l/r  as  a  function  of 
y/r  at  y/r  =  0.  The  alternative  assumption  introduced  by  Von  Kaiman 
regarding  the  mixing  length  term  is  based  on  considerations  of 
mechanical  similitude  of  the  turbulent  flow  pattern  from  one  point 
to  another.  This  leads  to  the  relation 

f  =  k(du/dy)/(d2u/dy2)  .  (3) 

When  either  of  the  values  of  I  are  introduced  in  Eq.  (l),  along  with 
the  constant  shear  stress  assumption,  the  result  in  integrated  form 
is  obtained  by  introducing  boundary  conditions  corresponding  to  the 
edge  of  a  laminar  sublayer  between  the  body  surface  anu  the  primary 
turbulent  portion  of  the  boundary  layer.  This  assumption  may  be 
expressed  in  the  form 


yu*/vw  =  5au*/vw  =  8  , 


(4) 
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in  which  is  the  kinematic  viscosity  at  the  wall,  u*  =  p  is  the 
so-called  friction  velocity,  and  the  entire  quantity  in  Eq.  (4)  is 
Prandtl's  friction- distance  parameter.  Since  yu*/vy  is  a  sort  of 
Reynolds  Number,  it  is  assumed  that  its  value  s  at  y  =  5  is  a 
critical  value  corresponding  to  the  transition  from  the  laminar  flow 
in  the  sublayer  to  the  turbulent  flow  outside  it.  Experimental 
evidence  indicates  that  s  -  11.6. 

2.1  Development  of  the  Law  of  the  Wall 

Before  carrying  out  the  integration  of  Eq.  (l),  it  is  convenient 
to  introduce  the  nondimens ional  ratios 

<p  =  u/u*  and  rj  =  yu*/vw  .  (5) 

The  integration  of  Eq.  (l)  with  either  Eq.  (2)  or  (3)  for  l  then 
gives  a  result  in  the  form 

<p  =  A  log1()Ti  +  B  ,  (6) 

where  the  constants  depend  on  the  values  assumed  for  k  and  s.  In  the 
case  of  an  incompressible  flow  without  heat  transfer,  the  integration 
yields  a  result  in  the  form 


9  -  -j^n  +  'n(v)j  .  (7) 

This  result  involves  the  assumption  that  the  velocity  distribution 
within  the  laminar  sublayer  is  linear,  that  is 

u  =  V/5s  »  (0) 
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being  the  so-called  wall  velocity  or  value  of  u  at  the  outer  edge 
of  the  sublayer.  The  shear  stress  within  the  sublayer  is  then  given 
as 


T 

O 


} 


(9) 


or,  after  eliminating  the  wall  velocity  and  replacing  /t Jq  by  u*, 
which  is  the  friction  velocity,  the  sublayer  velocity  distribution 
is  given  as 


u 

u* 


J2Z 

v 


(10) 


In  connection  with  the  turbulent  velocity  distribution 
represented  by  Eq.  (7),  the  values  of  the  constants  using  k  =  0.4 
and  s  =  11.6  are 


A  =  5.75  and  B  =  5.50  .  (ll) 

Schlichting  [Ref.  ( l) ]  suggests  a  slight  empirical  modification  of 
these  values  for  flat  plates  to 

A  =  5.85  and  B  =  5-56  ,  (12) 

00  that  Eq.  (6)  finally  takes  the  fora 

<P  =  5.85  log  tj  +  5.56  .  (15) 

This  result,  along  with  Eq.  (10)  for  the  sublayer  velocity 
distribution,  is  usually  referred  to  as  the  "law  of  the  wall." 

A  plot  of  cp  as  a  function  of  q  is  presented  in  Fig.  1  in 
semilogarlthnic  form,  based  on  Eqs.  (10)  and  (15). 
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When  compared  with  experimental  evidence,  two  deviations  from 

the  analysis  on  which  Fig.  1  is  based  are  immediately  evident.  The 

transition  from  the  laminar  sublayer  to  the  main  turbulent  flow  is 

not  Instantaneous  as  assumed  in  the  theory  at  tj  =  11.6,  but  rather 

takes  place  over  a  finite  portion  of  the  complete  boundary  layer. 

The  other  difference  is  found  in  the  outer  portion  of  the  boundary 

layer  where  the  values  of  <p  tend  to  be  somewhat  larger  than  those 

predicted  by  the  analysis.  A  good  comparison  with  experimental 

data  for  incompressible  flow  in  smooth  pipes  is  found  in  Fig.  20.4 

of  Ref.  (l).  Similar  plots  can  also  be  drawn  for  flat  plate  boundary 

layers,  although  in  this  case  a  complete  analysis  requires  a 

determination  of  local  shear  stress  t  in  order  to  evaluate  the 

o 

factor  u*.  Such  comparisons  will  be  given  later  in  connection  with 
compressible  flows. 

2.2  Determination  of  Skin  Friction  Coefficients 


The  results  contained  in  the  "law  of  the  wall"  equation  may  be 
used  as  a  basis  for  calculating  the  values  of  both  local  and  mean 
skin  friction  coefficients.  The  local  value  at  a  point  on  the  plate 
at  a  distance  x  from  the  leading  edge  is  given  by  the  relation 

Cf  =  2tc/PU1  '  (lJ0 

while  the  Reynolds  Number  based  on  the  same  distance  is 

R  =  pu-jx/u  ,  (15) 

the  assumption  here  being  that  the  boundary  layer  is  completely 
turbulent  beginning  at  the  leading  edge.  If  these  values  are  used 
to  detenaine  the  shear  stress  and  the  friction- distance  parameter, 
it  is  found  that  the  skin  friction  coefficient  and  the  Reynolds 
Number  are  related  by  the  expression 


9 


1/yf^  =  1.7  +  4.15  log  (cf  R) 


(16) 


The  numerical  coefficients  in  Eq.  (l6)  are  baaed  on  the  values  of 
A  and  B  represented  by  Eq.  (13). 

If  the  value  of  the  mean  skin  friction  coefficient  Is  desired, 
representing  the  drag  force  on  one  side  of  the  plate.  It  may  be 
written  In  the  form 


CF  =  X  .  (17) 

The  determination  of  the  value  of  the  drag  force  D„  is  made  by 

it  r 

using  the  Von  Karman-Fohlhausen  momentum  Integral  formula.  For 
one  side  of  the  flat  plate  In  Incompressible  flow,  this  relation 
Is  simply 


D 


f 


utuj-ujay 


(18) 


Introduction  of  the  law  of  the  wall  velocity  distribution  for  u, 
carrying  out  the  integration  across  the  boundary  layer  at  a  particular 
value  of  x,  and  the  omission  of  negligible  terms  in  the  result  finally 
lead  to  the  so-called  Von  Karman- Schoenherr  formula,  which  is 

0.242/yc^,  =  log(CF  R)  .  (19) 


Equation  (19)  has  the  disadvantage  that  it  is  not  possible  to 
solve  it  explicitly  for  CF  as  a  function  of  R,  although  the  converse 
relation  is  readily  obtainable.  Schlichting  [Ref.  (l),  p.  602]  has 
therefore  proposed  an  alternate  relation  which  closely  approximates 
Eq.  (19)  over  a  wide  range  of  Reynolds  Numbers  and  is  of  the  form 
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CF  =  0.455/(108  R) 


(20) 


2.58 


If  the  boundary  layer  has  an  initial  segment  of  laminar  flow, 
followed  by  a  transition  to  turbulent  flow,  an  analysis  by  Prandtl 
and  Schlichting  [Ref.  (l),  p.  602]  leads  to  a  relation  of  the  form 

Cp  =  0.455/{log  R)2*58  -  K/R  .  (21) 

This  result  is  based  on  the  combination  of  the  Blasius  formula  for 
a  laminar  incompressible  boundary  layer, 

Cp  =  1.328/(R)1/2  ,  (22) 

with  the  turbulent  skin  friction  relation  of  Eq.  (20).  The  value 
of  the  coefficient  K  in  Eq.  (2l)  depends  upon  the  Reynolds  Number 
at  which  the  transition  to  turbulent  flow  begins.  This  in  turn 
depends  on  the  shape  of  the  plate  leading  edge,  the  level  of  free 
stream  turbulence  in  the  case  of  a  wind  tunnel  test,  and  the  rough¬ 
ness  of  the  plate  surface,  particularly  in  the  vicinity  of  the 
leading  edge.  A  typical  value  of  this  critical  Reynolds  Number  is 
5  X  10 5,  which  leads  to  a  value  of  K  =  1700.  The  results 
represented  by  Eqs.  (20),  (2l),  and  (22)  are  shown  in  Fig.  2  in  a 
typical  logarithmic  plot. 

2.3  The  Velocity  Defect  Law 

As  already  mentioned,  the  law  of  the  wall  relation  of  Eq.  (6) 
does  not  agree  too  well  with  experimentally  determined  velocity 
distribution  data  for  large  values  of  n .  This  is  due  primarily 
to  the  fact  that  the  law  of  the  wall  was  derived  on  the  basis  of 
conditions  in  the  immediate  vicinity  of  the  wall  or  boundary  surface. 
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If  Instead  the  Integration  of  the  basic  relation  of  Eq.  (l)  is 
carried  out  using  conditions  at  the  outer  edge  of  the  boundary  layer, 
i.e.. 


y  =  6  ,  u  =  and  bn/by  =  0 

a  relation  of  the  form 


(23) 


<PL  -  1>  -  -A  log(y/6)  -  0  (24) 

is  obtained.  The  value  of  the  constant  A  is  found  to  be  5*75,  as 
in  the  case  of  the  law  of  the  wall,  and  this  value  holds  for  both 
pipe  and  boundary  layer  flows.  The  value  of  the  second  constant,  C, 
may  vary  widely  depending  on  the  nature  of  the  boundary  conditions. 

For  smooth  pipe  and  channel  flow,  a  value  of  C  *  -0.85  appears  to 
be  satisfactory,  while  for  boundary  layer  flows,  C  =  -3.25.  The 
relationship  of  Eq.  (24)  is  frequently  referred  to  as  the 
"velocity  defect  law,"  since  the  difference  -  qp  is  essentially 
the  difference  between  the  velocity  outside  the  boundary  layer  and 
the  value  within  it  for  values  of  y  approaching  the  boundary  layer 
thickness  6.  Both  values  are  of  course  expressed  as  ratios  to  the 
friction  velocity,  u*. 

In  comparing  the  law  of  the  wall  and  the  velocity  defect  law, 
it  should  be  noted  that  the  former  expresses  cp  as  a  function  of  the 
friction-distance  parameter  q  =  yu*/v,  while  the  latter  gives  <p  as 
a  function  of  y/6.  Thus,  if  both  laws  are  to  be  shown  graphically 
on  the  same  plot,  it  is  necessary  to  change  from  one  independent 
variable  to  the  other.  This  change  may  be  accomplished  by 
writing 


1 

6 


ZH!  JL  =  JJL) 

v  u*5  ™\u*5/ 
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The  friction  velocity  is  then  written  out 
the  shear  stress  is  expressed  in  terms  of 
coefficient,  c^,  with  the  result  that 

Thus,  it  is  necessary  to  have  information 
of  the  boundary  layer  thickness  6,  or  its 

ul® 

Rg  =  — — ,  and  also  the  value  of  c^,  in  order  to  accomplish  the 
desired  transfer.  In  Ref.  (2),  Schultz-Grunow  used  extensive  data 
from  measurements  of  the  turbulent  boundary  layer  on  a  flat  plate 
to  develop  this  idea  further,  including  expressions  for  R^  and  for 
cf .  He  then  plotted  his  velocity  profile  data  both  as  a  function 
of  T),  as  well  as  a  function  of  y/8.  His  plot  in  the  latter  form 
is  shown  in  Fig.  3  and  clearly  illustrates  the  significance  of  the 
velocity  defect  law  as  it  applies  to  the  outer  portion  of  the 
boundary  layer. 

2.4  The  Law  of  the  Wake 


as  u*  =  yj  t  Jp  and 
the  local  skin  friction 


(25) 

available  as  to  the  value 
Reynolds  Number 


An  extension  of  the  velocity  defect  concept  has  been  developed 
by  Coles  [Ref.  (3)],  a  brief  summary  of  which  is  given  in  an  article 
by  W.  C.  Reynolds  [Ref.  (4)].  The  results  of  Coles'  analysis  as 
summarized  in  Reynolds'  article  are  given  in  the  so-called 
"law  of  the  wake"  which  may  be  written  in  the  form 


<f>l  -  9  =  -  ^  *n(y/5)  +  £  [2-w(y/&)]  .  (26) 


The  coefficient  n  is  described  as  a  "wake  factor"  which  is  related 
to  the  shape  parameter  H^g  and  the  local  skin  friction  coefficient  cf. 


15 


The  shape  parameter  used  here  Is  the  usual  value  detemlned  by  the 
ratio  of  the  displacement  and  manentun  thicknesses,  that  Is, 

\2  =  V62  *  (27) 

The  value  of  k  Is  the  usual  Prandtl-Nikuradse  mixing  length  factor, 
taken  as  0.4  In  the  present  discussion.  Finally  the  function  w(y/6) 

Is  Coles'  wake  function  which  he  originally  gave  in  numerical  form 
and  which  has  been  approximated  by  various  empirical  relations, 
typical  of  the  latter  are  the  transcendental  relation 

w(y/8)  *  1  -  cos  (jtrj)  (28) 

used  by  Hlnze  and  Spalding  and  the  polynomial  form  suggested  by 
Rotta  [Ref.  (5)], 

v(y/B)  =  39t)3  -  I25i4  +  I85ri5  -  133tj6  +  W  .  (29) 

Rotta  indicates  that  the  law  of  the  wake  may  be  used  to  obtain  a 
relation  for  the  local  skin  friction  coefficient  In  the  form 

=  i  inpT) +  c'  +  K'  >  (50) 

where  5^  is  the  displacement  thickness  of  the  boundary  layer.  He 
Indicates  that  C'  is  the  usual  constant  value  of  3*2  appearing  In 
the  law  of  the  wall  relation,  while  the  value  of  the  second  constant 
Is  Indicated  to  be  obtainable  as  a  function  of  the  wake  factor  n. 

The  numerical  values  of  Coles'  laws  were  given  in  Ref.  (3)  for 
the  wall  region  as  well  as  for  the  outer  portion  of  the  boundary 
layer.  These  values  were  based  on  an  extensive  analysis  of  all 
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experimental  data  then  available  and  resulted  In  the  values  given  here 
In  Table  I.  The  basic  relations  to  which  these  data  apply  are  the 
law  of  the  wall,  written  in  the  form 

<P  «  fjh)  ,  (31) 

and  the  velocity  defect  law, 

<PX  -  9  =  f2(^)  •  (32) 

It  is  now  quite  generally  accepted  practice  to  use  these  values,  or 
empirical  relations  based  on  them,  in  lieu  of  expressions  such  as 
Eqs .  (13)  and  (24)  which  were  derived  from  mixing  length  theory. 

In  connection  with  these  last  two  references,  it  should  be 
noted  that  the  paper  by  Reynolds  appears  in  the  Proceedings  of  an 
invited  conference  held  at  Stanford  University  in  the  summer  of  1968, 
in  which  an  international  group  of  workers  in  the  turbulent  boundary 
layer  field  who  had  developed  various  calculation  methods  presented 
their  procedures  as  applied  to  a  number  of  selected  incompressible 
flow  problems.  These  methods  were  subjected  to  a  close  critical 
evaluation  of  their  effectiveness,  so  that  this  publication  represents 
a  comprehensive  summary  of  the  state-of-the-art  at  the  time  of  its 
publication.  The  Proceedings  also  includes  a  paper  by  Rotta  which 
presumably  covers  some  of  his  later  work  in  this  field  beyond  that 
represented  by  Ref.  (5).  The  methods  considered  by  the  Stanford 
conference  were  limited  to  incompressible  flows,  but  included  some 
considerations  of  roughness  and  external  pressure  gradient  effects. 
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Ill 


TURBULENT  BOUNDARY  LAYERS  IN  COMPRESSIBLE  FLOW 


The  development  of  interest  in  the  boundary  layer  problems  of 
compressible  fluids  is  actually  one  of  long  standing,  as  evidenced 
by  the  application  to  such  flows  in  high-speed  fluid  machinery 
involving  gaseous  fluids  and  in  aeroballistics .  Active  research  in 
this  area  did  not  progress  substantially  until  the  development  of 
high-speed  aircraft  and  guided  missiles;  in  fact,  it  was  in  connection 
with  design  problems  In  the  latter  field  that  the  ARL  boundary  layer 
research  program  was  Initiated  late  in  1947  •  A  dual  approach  was 
formulated,  involving  an  extension  of  incompressible  flow  boundary 
layer  theory,  as  well  as  an  experimental  program  of  velocity  distri¬ 
bution  measurements  and  skin  friction  determinations.  The  next 
sections  of  this  report  will  present  a  summary  of  same  of  the  more 
significant  work  concerned  with  the  analytical  phase  of  the  program. 

3.1  Boundary  Layer  Theory  for  Compressible  Fluids 

An  early  approximate  method  of  estimating  skin  friction  in  a 
compressible  fluid  was  devised  by  Von  Karman  [Ref.  (6)],  based  on 
the  concept  of  replacing  the  free  stream  values  of  viscosity  and 
density  by  values  corresponding  to  the  wall  temperature.  The 
temperature  itself  is  determined  by  writing  the  energy  equation 
without  dissipation  for  a  local  point  in  the  flow  and  the  free 
stream,  that  is, 


C  T 
P 


°PTl  +  T 


C  T 
P 


(55) 


Preceding  page  blank 
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The  Introduction  of  the  free  stream  Mach  Number  makes  it  possible 
to  solve  for  the  ratio  of  local  and  wall  temperatures,  the  result 
being 


When  applied  to  free  stream  conditions,  u  =  and  T  =  T1, 
Eq.  (54)  becomes 


(35) 


The  introduction  of  this  value  in  Eq.  (34)  finally  gives  for  the 
ratio  of  local  and  vail  temperatures  the  expression 


where 


a  - 


m 

1  +  m 


(37) 


and 


m  =  **1  *  (38) 

All  of  this  analysis  is  based  on  the  assumption  that  there  is  no 
significant  heat  transfer  at  the  wall  and  no  dissipation  of  energy 
within  the  boundary  layer. 
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The  density  is  next  assumed,  to  be  inversely  proportional  to  the 
temperature,  while  the  viscosity  relation  may  be  approximated  by  a 
simple  exponential  relation  of  the  form 


±L  = 


(39) 


with  co  =  0.768  for  air.  The  introduction  of  the  wall  values  in  the 
simple  boundary  layer  momentum  theory,  based  on  a  one- seventh  power 
velocity  distribution,  yields  the  result  that 


(40) 


The  same  procedure  can  also  be  carried  out  on  the  basis  of  the 
law  of  the  wall  skin  friction  formula,  Eqs.  (6)  and  (19),  with  the 
result  that 


where  Tw/P1  is  calculated  from  Eq.  (35).  Either  Eq.  (40)  or  (4l) 
may  be  used  to  calculate  the  ratio  of  the  compressible  and 
incompressible  values  of  the  skin  friction  coefficients  as  functions 
of  free  stream  Mach  Number,  the  results  being  shown  in  Fig.  4. 

A  more  rigorous  analysis  of  the  compressible  turbulent  boundary 
layer  was  carried  out  by  Wilson  [Ref.  (7)]  as  the  initial  analytical 
part  of  the  ARL  program.  He  developed  his  analysis  by  beginning  with 
the  Prandtl  mixing  length  relation  of  Eq.  (l)  and  the  Von  Kaman 
mixing  length  formulation  of  Eq.  (3),  with  the  density  and  viscosity 
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terns  being  expressed  in  terms  of  the  wall  temperature  as  in  the 
preceding  case,  but  with  the  modification  of  introducing  a 
temperature  recovery  factor, 


r 


so  that 


T[  -  1  +  r(¥)  A  <«) 

replaces  the  relation  of  Eq.  (35).  The  integration,  which  is  rather 
lengthy,  finally  leads  to  a  law  of  the  wall  for  a  compressible 
boundary  layer  in  the  form 


The  values  of  <p  and  cp^  are  the  velocity  ratios,  u/u*,  within  and  at 
the  outer  edge  of  the  boundary  layer,  as  in  the  incompressible  case, 
while  8  again  is  the  value  corresponding  to  the  outer  edge  of  the 
laminar  sublayer.  When  converted  to  logarithmic  form,  Eq.  (44) 
becomes 
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but  instead  of  evaluating  the  right-hand  side  of  this  expression  in 
terms  of  tj,  present  practice  is  to  assume  that  the  right  side  may  be 
regarded  as  a  function  of  tj  which  is  best  defined  by  the  numerical 
data  established  by  Coles  [Ref.  (3)],  which  have  been  used  as  a 
basis  for  plotting  the  curve  shown  in  Fig.  4.  It  is  significant  to 
note  that  when  plotted  in  this  form,  the  curve  of  Fig.  4  is 
independent  of  Mach  Number.  Also  the  transition  from  the  laminar 
sublayer  to  the  main  turbulent  portion  of  the  boundary  layer  is 
included  in  Coles'  data.  Included  in  Fig.  4  are  some  experimental 
measurements  made  at  ARL  for  several  Mach  Numbers  and  for  slight 

variations  in  the  momentum  thickness  Reynolds  Numbers,  R.  . 

2 

These  data  are  taken  from  the  work  of  Fenter  and  Stalmach  [Ref.  (8)] 
and  show  the  trend  away  from  the  law  of  the  wall  into  the  law  of  the 
wake  for  the  larger  values  of  t). 

Wilson  did  not  undertake  a  comparison  with  either  the  law  of 
the  wall  or  the  law  of  the  wake  using  his  velocity  profile  measurements. 
He  did,  however,  transform  the  results  represented  by  Eq.  (46)  into  a 
relation  between  mean  skin  friction  coefficient  and  Reynolds  Number. 

This  was  accomplished  by  first  introducing  values  of  C_  and  R„,  both 
based  on  density  and  viscosity  values  at  the  plate  surface.  A  further 
transformation,  introducing  free  stream  values  of  these  latter 
quantities  and  integrating  along  the  plate  by  means  of  the  momentum 
equation,  finally  led  to  the  result  that 


0.242 >/(l+rm)Cy  =  log  (C^)  -  <0  log  (l+im) 


(47) 


Attention  is  invited  to  the  fact  that  0  and  m  are  the  Mach  Number 
functions  defined  by  Eqs.  (37)  and  {38).  It  should  also  be  noted 
that  Eq.  (47)  is  based  on  Wilson's  original  form  as  given  by  Eq.  (46), 
rather  than  on  Coles'  data  for  f(T]). 
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The  result  of  Eq.  (1*7)  closely  resembles  the  Von  Kaiman- Schoenerr 
formula  of  Eq.  (19),  except  for  the  modifications  represented  by  the 
presence  of  the  Mach  Nunber  functions,  o  and  m,  and  the  temperature 
recovery  factor  r.  An  important  feature  of  these  results  may  be 
demonstrated  by  a  rearrangement  of  Eq.  (47)  into  the  form 


_ 0.242 _ 

: 

[«(«m)/(»iii‘1oVa)  ] 


=  log 


aj  1+m) 


rain  g 
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1  1/2' 


(l+m) 


(l+o>) 


or  by  replacing  the  two  Mach  Number  functions  by 


(48) 


X  = 


(l+m)Cp 
(sin_1o1/2) 


(^9) 


and 


(50) 


Eq.  (48)  may  be  written  in  the  condensed  form 


0.242  _ 

~p72  ‘  106  ** 


(51) 


The  important  feature  of  this  result  is  that  a  plot  of  X  as  a 
function  of  Y  is  obtained  as  a  single  universal  curve  which  is  the 
same  for  all  Mach  Numbers.  For  numerical  calculations,  it  is 
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necessary  to  solve  for  Y  as  a  function  of  X,  or  in  effect, 
Reynolds  Number  as  a  function  of  skin  friction  coefficient.  This 
form  is 


log  Y  =  "  log  X  .  (52) 

A  further  simplification  in  the  calculation  procedure  is  obtained 
by  separating  out  the  portions  of  X  and  Y  which  are  dependent  on 
Mach  Number  only.  This  is  accomplished  by  writing 

X  =  CF  X'  (55) 

and 

*  -  V  .  (54) 

where 


Y  = 


stn-V/2) 

o(l+rm)^+a3 


(56) 


Graphs  showing  the  variation  in  the  several  Mach  Number  functions  are 
presented  in  Figs.  5  and  6.  In  the  first  of  these  figures  the  value 
of  the  term  1  +  rm  is  shown  for  two  values  of  the  recovery  factor, 
r  =  1.00  and  r  =  0.88,  in  order  to  illustrate  the  small  effect  of  this 
factor.  In  Fig.  6  are  shown  the  values  of  l/(l  +  m)  and  (l  +  m)1+aJ 
for  a>  =  O.768,  the  additional  Mach  Number  factors  required. 
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Figures  7  and  8  show  the  values  of  the  parameters  X'  and  Y',  as  well 
as  their  reciprocals,  all  as  functions  of  Mach  Number.  Making  use 
of  all  of  these  factors,  it  is  now  possible  to  compute  the  data 
represented  by  the  curve  shown  in  Fig.  9*  This  is  the  uni  versed 
relation  between  skin  friction  coefficient  and  Reynolds  Number 
described  by  Eqs.  (5l)  or  (52).  The  curve  of  Fig.  9  was  used  as  a 
basis  for  determining  the  actual  values  of  Cj,  as  a  function  of 
Reynolds  Number,  R^,  for  three  different  values  of  the  free  stream 
Mach  Number,  that  is,  =  0  corresponding  to  incompressible  flow 
and  two  supersonic  values,  =  5  and  =  10,  plotted  in  Fig.  10. 
Finally  curves  are  plotted  in  Fig.  11  showing  the  ratio  of  the 
compressible  and  incompressible  values  of  Cj.,  this  being  done  for 
an  assuned  Reynolds  Number  of  10^.  Similar  curves  are  also  shown 
based  on  Von  Harman' s  early  analysis  which  led  to  Eq.  (41),  the 
two  values  of  the  recovery  factor  again  being  considered.  All  of 
these  curves  are  shown  as  functions  of  Mach  Number  in  Fig.  11. 

In  this  connection  reference  should  also  be  made  to  a  more 
recent  work  by  Wilson  [Ref.  (8)]  which  includes  seme  of  the  previous 
material  on  smooth,  adiabatic  plates  but  also  adds  more  recent 
work  on  heat  transfer  and  roughness  effects. 
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IV.  THE  EFFECTS  OF  SURFACE  ROUGHNESS 

The  determination  of  the  effects  of  roughness  on  the  character 
of  the  boundary  layer  flow  is  a  problem  of  long  standing.  One  of 
the  earliest  systematic  investigations  was  the  work  of  Nikuradse 
which  actually  involved  pipe  flow  rather  than  boundary  layer  flow. 

In  both  cases,  however,  roughness  patterns  may  be  classified  as 
(l)  unifonnly  distributed,  (2)  distribution  of  isolated  roughness 
elements,  and  (3)  various  special  patterns  which  do  not  fall  into 
either  of  the  first  two  categories.  Nikuradse 's  work  on  pipe  flow 
was  accomplished  by  coating  the  interior  wall  with  sand  grains  of  a 
specified  average  size  and  then  determining  the  longitudinal  pressure 
drop  and  the  velocity  distribution  across  the  pipe.  This  variety  of 
roughness  falls  in  the  category  of  the  uniformly  distributed  type, 
although  the  precise  shape  of  the  roughness  particles  could  have 
varied  substantially.  In  seme  of  the  early  work  done  at  ARL  on  the 
boundary  layer  roughness  problem,  uniformly  distributed  roughness 
was  described  as  being  of  the  "Rocky  Mountain  type,"  exemplified  by 
the  distributed  sand  grains,  or  of  the  "Mole  Hill  type"  when  the 
roughness  particles  are  smooth  and  regular,  such  as  would  be  the 
case  when  small  spherical  beads  constituted  the  roughness  elements. 

4.1  Roughness  Effects  in  Incompressible  Flow  in  Pipes 

A  comprehensive  sunnary  of  Nikuradse' s  work  with  what  was 
essentially  incompressible  flow  is  given  in  Ref.  (l),  Chaps.  XX- f 
and  XXI- c.  The  mechanism  of  thfe  flow  has  been  well  demonstrated 
in  terms  of  boundary  layer  concepts  by  recognizing  the  fact  that  for 
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a  given  roughness  particle  size  and  considering  flows  with  increasing 
Reynolds  Number,  the  boundary  layer  initially  completely  covers  the 
roughness  projections  so  that  the  friction  factor  as  a  function 
of  Reynolds  Number  follows  closely  the  variation  found  for  smooth 
surfaces.  As  the  Reynolds  Number  Increases,  there  Is  a  steady 
decrease  In  boundary  layer  thickness  until  a  point  is  reached  where 
the  roughness  elements  extend  through  It  Into  the  main  turbulent 
portion  of  the  flow.  There  is  thus  an  augmentation  of  the  turbulent 
flow  losses  and  a  gradual  leveling  off  of  the  friction  factor  for 
further  Increases  in  Reynolds  Number.  The  conditions  at  which  this 
situation  develops  have  been  referred  to  as  “admissible"  or 
"threshold  roughness."  There  Is  also  a  change  In  the  velocity 
distribution  across  the  pipe,  characterized  by  a  relation  of  the  fora 


u 

u* 


=  5.75  log 


+  B 


(57) 


where 


(1) 

u  Is 

(2) 

u#  is 

(5) 

k  is 
s 

(M 

B  is 

the  pipe  wall, 


Number,  u*k  /v.  The  nature  of  the  variation  In  B  is 

described  by  the  plot  of  Pig.  12. 

On  the  basis  of  these  data,  smooth  pipe  flow  may  be  considered  as 

corresponding  to  roughness  Reynolds  Numbers,  u*k  /v<5,  in  which  case 

s 

B  is  essentially  a  linear  function  of  u*k  /v  of  the  fora 

s 


B  =  5.56  +  5.85  log  (u*k  /v) 

s 


(58) 
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There  next  follows  a  transition  zone  in  which  B  rises  to  maximum  and 
then  decreases  slightly  to  about  8.5  in  the  completely  rough  zone. 
These  three  zones  are  then  characterized  by  the  following  values  of 
the  roughness  Reynolds  Number: 

Hydraulically  smooth:  (u*k  /v)<5 

8 

Transition :  5<(  u*k  /v  )<70 

s 

Completely  rough:  (u*k  /v)>70 

s 

Additional,  studies  of  rough  pipes  have  shown  that  in  the  case  of 
Isolated  roughness  patterns  with  roughness  elements  located  at 
discrete  points,  and  for  commercial  pipe  of  various  types,  the 
roughness  can  be  expressed  in  terns  of  an  equivalent  sand  grain 
roughness;  these  results  are  alf;o  summarized  in  Ref.  (l). 

k.2  Roughness  Effects  in  Incompressible  Flow  on  Flat  Plates 

Prandtl  and  Schlichting  [Ref.  (lO)]  made  a  very  useful 
extension  of  Nlkuradse's  pipe  flow  research  by  applying  his  results 
to  flat  plates  as  affected  by  sand  grain  roughness.  They  also 
demonstrated  that  the  concept  of  admissible  roughness  continued  to 
be  applicable,  with  the  value  of  the  admissible  sand  grain 
dimension  given  by 


s  100v^  •  (59) 

In  addition  to  an  elaborate  set  of  charts  for  estimating  roughness 
effects  on  flat  plate  skin  friction,  interpolation  formulas  were 
suggested  for  local  and  mean  skin  friction  coefficients.  These 
relations  are 
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(60) 


c 


f 


2.87  +  1.58  log 


-2.5 


and 


C 


P 


1.89  +  1.62 


(61) 


In  which  l  Is  the  total  length  of  the  plate  and  x  Is  the  distance 
along  the  plate  from  Its  leading  edge  to  the  point  where  the  local 
skin  friction  coefficient  is  to  be  calculated.  These  relations  are 

o  f. 

considered  valid  for  10  <  l /is.  <  10. 

4.3  Roughness  Effects  in  Compressible  Plow  on  Plat  Plates 

In  a  rather  detailed  study  of  the  effects  of  roughness  on  flat 
plates,  Fenter  [Ref.  (ll) ]  suggested  that  the  factor  B  in  Eq.  (57) 
could  be  conveniently  approximated  by  a  series  of  linear  semilogarithmic 
relations.  Using  the  type  of  roughness  function  for  flat  plates,  the 
law  of  the  wall  is  developed  in  terms  of  the  mean  roughness  height  kr, 
and  the  corresponding  value  of  the  friction-distance  parameter  is  then 
written  as 


\  -  kruVvw  •  (62) 

The  law  of  the  wall  for  a  roughened  plate  may  then  be  shown  to  take 
the  form 


cp  =  Ar  log(y/kr)  +  Br 


(63) 


28 


The  first  coefficient,  A^,  is  equal  to  A  as  in  Eq.  (6)  for  smooth 
surfaces;  that  is,  Af  =  A  =  5.75.  Equation  (63)  may  be  expanded 
into  the  form 

cp  =  Ar  log  y  -  Ar  log  kr  +  Br  (64) 

in  which  the  roughness  effect  is  now  represented  by  the  last  two 
terms;  that  is 

b;  -  Br  -  Ar  log  tr  .  (65) 

In  order  to  Include  the  smooth  plate  case  in  this  analysis,  Fenter 
replaced  Eq.  (65)  by  the  expression 

B^  =  B  -  f(i]r)  ,  (66) 

where  B  =  5.50  is  the  value  previously  used  in  the  incompressible 
smooth  surface  case  and  f(rjr)  is  a  new  roughness  function  which,  on 
the  basis  of  the  linear  approximations  referred  to  previously,  may  be 
written  in  the  form 

f( Tlr)  =  a  log  rjr  -  b  .  (67) 

This  modified  roughness  function  was  assumed  to  be  determined  directly 
from  Nikuradse's  incompressible  pipe  flew  data  given  in  Fig.  12.  The 
values  of  f(T)r)  for  the  three  flew  regimes,  along  with  the  corresponding 
forms  of  Eq.  (67), are  plotted  as  functions  of  log  tj  in  Fig.  13.  The 
values  of  the  coefficients  a  and  b  and  the  ranges  of  values  of  tj  are 
given  in  Table  I.  It  should  be  noted  that  the  value  of  tj  denoting  the 
demarcation  between  the  transition  and  the  fully  rough  regimes  has 
been  taken  as  tj^  =  100,  instead  of  the  slightly  lower  value  of  70 
suggested  by  Schlichting.  A  similar  treatment  of  the  roughness  function 
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was  also  used  by  Young  [Ref.  (12)}  in  a  later  phase  of  the  AKL  program. 
Involving  heat  transfer  as  well  as  roughness  but  subdividing  the 
roughness  zones  Into  as  many  as  five  segments,  as  shown  In  Fig.  12. 

It  Is  probable  that  a  single  functional  relation  could  be 
established  for  f(  r)r)  covering  the  full  range  of  tj  values,  but 
this  has  not  been  undertaken  up  to  the  present  time. 


less  and  Heat  Transfer  Effects  In  Compressible  flow 


Fenter  Included  In  his  work  presented  In  Ref.  (ll)  an  extension 
of  the  theoretical  analysis  leading  to  integrated  forms  of  the 
velocity  profile  results  which  finally  yield  relations  for  the  local 
and  mean  skin  friction  coefficients  as  functions  of  Reynolds  Number. 
Actually  the  theory  was  developed  In  a  rather  general  form  so  as 
to  Include  simultaneously  the  effects  of  both  roughness  and  heat 
transfer.  While  the  basic  approach  was  comparable  to  that  used  by 
Wilson  [Ref.  (7)]  in  his  earlier  treatment  of  ccmpresslble  flow  on 
a  smooth  adiabatic  plate,  the  details  are  considerably  more  formidable. 
After  a  fairly  extensive  analysis,  Fenter' s  results  indicate  that  the 
law  of  the  wall  velocity  profile  may  be  written  as 


for  a  fully  rough  surface.  The  factors  cp*,  0,  and  C  appearing  In 
Eq.  (68)  are  defined  as 
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(70) 


fo  (2ocp/cp1-X) 
kcpx  A2+i| -ff 


(71) 


The  value  of  or  is  the  Mach  Number  function  previously  introduced  in 
Eg.  (36);  that  is. 


while 


(72) 


(73) 


The  corresponding  fom  of  the  velocity  defect  law  is 

<PJ  -  <P*  =  F(y/6)  ,  (74) 


where  F(y/b)  is  based  on  the  Coles'  function  given  in  numerical 
fom  in  Table  II. 


A  similar  approach  to  the  boundary  layer  problem  for  rough 
surfaces  with  heat  transfer  was  developed  by  Van  Driest  [Ref.  (13)] 
using  Prandtl's  mixing  length  theory  with  Von  Karman’s  relation  for 
the  mixing  length  tern  itself.  These  results  were  utilized  by  Young 
[Ref.  (12)]  to  provide  a  basis  of  comparison  between  theory  and 
experiment  in  connection  with  an  experimental  study  of  Veegroove  type 
roughness  coupled  with  heat  transfer  on  a  flat  plate.  Van  Driest' s 
method  of  calculation  was  also  included  in  an  extensive  comparison  of 
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theories  with  experimental  data  by  Spalding  and  Chi  [Ref.  (l4)],  the 
conclusion  being  that  Van  Driest  *s  approach  appeared  to  be  the  most 
satisfactory.  The  only  Important  restriction  in  the  development  of 
Van  Driest* s  method  was  that  the  Prandtl  Number  was  assumed  to  be 
constant  and  equal  to  unity. 

By  making  use  of  the  momentum  integral  relation,  both  Renter's 
and  Van  Driest *s  results  may  be  used  to  derive  an  expression  for  the 
local  skin  friction  coefficient  in  terms  of  the  momentum  thickness 
Reynolds  Number.  In  the  form  given  by  Penter,  this  result  is 


nr  t 

V  TT  “7=  =  log  RS  -  -F  +  log  ST  +  2.90 

v/rt  V  tj  . .  n  °0  Jo 


The  skin  friction  coefficient  in  Eq.  (75)  is  based  entirely  on  free 
stream  conditions,  that  is 


Similarly,  the  momentum  thickness  Reynolds  Number  is 


P1  U1  S2 


The  other  factors  in  Eq.  (75)  are 


(76) 


(77) 


52 


which  on  comparison  with  Eq.  (69)  for  c p*  will  be  seen  to  be 
equivalent  to  cp*  >/o/cp^  when  cp  =  qp^.  The  factor  f  in  Eq.  (75)  is 
the  roughness  function  f(nr)  previously  discussed  and  represented 
by  the  plot  of  Fig.  13 .  Finally  it  should  be  noted  that  the  value 
of  tj  may  be  determined  from  the  relation 


Rr  being  the  Reynolds  Number  based  on  the  mean  height  of  the 
roughness  elements  and  the  free  stream  velocity  and  viscosity.  The 
numerical  factors  in  Eq.  (75)  have  been  modified  slightly  from  those 
calculated  from  the  law  of  the  wall,  to  obtain  somewhat  improved 
agreement  with  experimental  results. 

For  engineering  purposes,  it  is  convenient  to  have  available 
plots  of  mean  skin  friction  coefficients  as  functions  of  the  normal 
Reynolds  Number  based  on  distance  along  the  plate  surface.  Such 
plots  were  first  prepared  by  Moody  [Ref.  (15)]  for  roughened  pipes, 
with  the  roughness  Reynolds  Number  based  on  the  equivalent  send 
grain  roughness  dimension  as  a  parameter.  Similar  plots  for 
roughened  flat  plates  are  given  by  Prandtl  and  Schlichting 
[Ref.  (10)]  based  on  the  assumption  that  the  flow  is  turbulent 
from  the  leading  edge  rearward.  An  additional  presentation  of  flat 
plate  data  for  a  range  of  Mach  Numbers  up  to  5.0  has  been  developed 
by  Clutter  [Ref.  (l6)]  and  is  also  presented  by  Wilson  [Ref.  (9)]. 
Clutter's  calculations  are  based  on  Van  Driest' s  analysis  for  the 
smooth  plate  and  on  a  simplified  treatment  of  the  roughness  problem 
developed  by  Liejmann  and  Goddard  [Ref.  (17)]. 

The  inclusion  of  heat  transfer  effects  on  boundary  layer  flows 
has  already  been  discussed  to  some  extent  in  connection  with  the 
preceding  analysis.  Only  the  simplest  case  of  a  variation  in  plate 
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or  wall  temperature  as  compared  to  the  adiabatic  or  zero  heat  transfer 
case  will  be  considered.  In  general  the  case  of  heat  transfer  to  or 
from  the  plate  may  be  described  In  terms  of  the  ratio  of  wall 
temperature  to  free  stream  temperature,  this  being  represented  by 
the  factor  T which  appears  in  the  equations  Just  presented.  No 
attempt  will  be  made  here  to  discuss  the  more  complicated  case  of  an 
ablating  surface  which  may  introduce  a  film  of  molten  plate  material 
along  with  the  possibility  of  chemical  reactions  taking  place  in  this 
region. 

In  addition  to  the  determination  of  local  and  mean  skin  friction 
coefficients  along  with  boundary  layer  velocity  profiles,  it  is 
usually  important  to  provide  a  means  for  the  calculation  of  heat 
transfer  rates.  This  evaluation  may  be  accomplished  in  terms  of  a 
variety  of  parameters  in  addition  to  the  flow  factors  of  Reynolds 
Number  and  Mach  Number.  The  latter  terms  are  primarily  related  to 
similarity  conditions  based  on  the  moment  un  or  Navier- Stokes 
equations,  while  the  heat  transfer  parameters  are  obtained  by 
similar  analyses  of  the  energy  equation.  The  first  of  these  parameters, 
the  so-called  Prandtl  Number,  is  defined  as 

P  =  HCp/x  .  (79) 

It  is  interesting  to  note  that  the  Prandtl  Number  is  dependent  only 
on  the  physical  properties  of  the  fluid  and  not  on  any  of  the  flow 
characteristics.  It  may  be  interpreted  physically  as  a  representation 
of  the  ratio  of  the  rate  of  diffusion  of  vorticity,  p/p,  to  the  rate 
of  diffusion  of  heat,  X/pc^,  where  X  is  the  thermal  conductivity  and 
^  is  the  specific  heat  at  constant  pressure.  The  thermal 
conductivity  factor  is  first  encountered  in  writing  down  the  basic 
Fourier  relation  for  convective  heat  transfer,  that  is. 
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where  q  is  the  heat  transfer  rate  per  unit  time  and  unit  area  of 
the  surface. 


Another  factor  of  importance  in  determining  heat  transfer  is 
the  combination  of  temperatures  known  as  the  Eckert  Number,  defined  as 


E  = 


2(VTl) 

T  -T, 
w  1 


(81) 


where  the  new  term,  T  ,  is  the  temperature  at  a  stagnation  point. 

s 

Assuming  fully  adiabatic  conditions,  the  value  of  T  is  determined  by 

s 

the  adiabatic  temperature  rise, 


T 

s 


(82) 


Equations  (8l)  and  (82)  may  be  combined  to  yield  a  relation  between 
the  Eckert  and  Mach  Numbers  of  the  form 


E  = 


2m 

(VV'i 


(83) 


The  heat  transfer  rate  determined  by  Eq.  (80)  may  alco  be 
expressed  in  terms  of  a  heat  transfer  coefficient  and  the  difference 
between  the  recovery  temperature  and  the  actual  wall  temperature. 

Thus 

»  *  »(Tr-Tu)  ,  (84) 

where  T^  corresponds  to  the  condition  that  =  0.  The 

combination  of  Eqs.  (80)  and  (84)  gives 
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(85) 


h(T  -T 
'  r  w 


) 


which  may  be  nondimenslonalized  by  writing  y  =  y 1 i,  with  I  as  a 

geometric  scale  factor,  and  T  =  T"(T  -T  )  for  the  local  temperature. 

r  w 

Equation  (85)  then  reduces  to  the  fonn 


(86) 


which  defines  a  new  nondimens lonal  heat  transfer  parameter  known 
as  the  Nusselt  Number.  Still  another  parameter  may  be  Introduced 
by  arbitrarily  inserting  the  Prandtl  Number  and  the  Reynolds  Number 
in  Eq.  (86).  This  procedure  shows  that  the  Nusselt  Number  may  then 
be  expressed  as 

(87) 


where  the  new  combination 


u. 


=  S 


is  known  as  the  Stanton  Number. 


(88) 


There  are  additional  heat  transfer  parameters  known  as  the 
Peclet  Number  and  the  Grashof  Number,  but  these  are  not  pertinent 
to  the  present  discussion.  A  comprehensive  summary  of  these  elements 
of  heat  transfer  theory  can  be  found  in  Ref.  (l),  Chap.  XII,  based 
on  laminar  boundary  layer  theory.  In  tne  case  of  turbulent  boundary 
layers,  the  same  factors  are  utilized,  its  being  necessary  to  make 
use  of  equivalent  or  eddy  viscosity  and  conductivity. 
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Approximate  solutions  to  an  important  group  of  gas  flow  problems 
may  be  obtained  by  assuming  that  the  Prandtl  Number  is  equal  to 
unity.  This  assumption  leads  to  considerable  simplification  of  the 
basic  equations  and,  in  the  case  of  air,  is  not  too  far  from  reality. 
For  air  over  a  moderate  range  of  temperatures,  the  Prandtl  Number 
has  a  value  between  0.71  and  0.72.  If  the  Prandtl  Number  is  assumed 
to  be  unity,  it  may  be  shown  that  the  temperature  within  a  boundary 
layer  is  given  by  the  Crocco  relation  as  a  function  of  local  velocity; 
that  is. 


TAX  =  TwAi  -  (W1)  Vu!  +  m(u/u1)(l-u/u1)  .  (89) 

For  the  more  general  case  where  the  Prandtl  Number  is  a  variable  and 
not  equal  to  unity,  Harkness  [Ref.  (l8)]  has  suggested  the  replace¬ 
ment  of  the  quadratic  relation  of  Eq.  (89)  by  a  cubic  of  the  form 

TAi .  +  (c+i)[(TrA1)  -  (y*i)] 

+  ^3.627(l-T +  2.62J  m  -  cj(T -  (T )][u/Ul)2  (90) 

-  2. 627^1  +  m  -  (Tr/T1)J(u/u:i)5 

The  coefficient  c  in  this  expression  is  a  function  of  the  Prandtl 
Number  of  the  form 


c 


mfP  -P 
w  w. 
ins. 

BW- . WJ 


(91) 
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where  P  and  P .  sure  respectively  the  Prandtl  Numbers  for  the 
w  ins 

actual  wall  and  for  Insulated  wall  conditions.  In  attempting  to 
arrive  at  a  skin  friction  evaluation,  Harkness  introduced  the 
assumption  that  the  laminar  sublayer  thickness  was  a  function  of 
the  heat  transfer  rate.  Based  on  a  limited  amount  of  experimental 
data,  this  relation  is 

s  =  sQ  +  6.6t  ,  (92) 

where  s  is  the  sublayer  parameter  6su*/v,  sq  a  u.6  is  its  value  for 
the  zero  heat  transfer  case,  and  t  is  a  temperature  factor  of  the 
form 

t  -  1  -  V*r  •  (93) 

The  integration  of  the  mixing  length  equation  using  Harkness 's 
modifications  leads  to  an  extremely  complicated  and  awkward  form. 

An  alternate  approach  using  a  simpler  quadratic  temperature  variation 
was  developed  by  Moore  [Ref.  (19)]  as  a  later  phase  of  the  ARL  program. 
The  relation  between  temperature  and  velocity  was  taken  in  the  form 


w 
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In  addition  Moore  also  used  Harkness's  relation  for  the  sublayer 
thickness  as  given  by  Eq.  (93)  f  with  the  result  that  the  law  of  the 
wall  relation  becomes 


Gt)  8in"1  (*•?) ' 6,6t  +  (s) fn (v)  =  tM  ■  ( 

The  right  side  of  this  expression,  f ( tj),  is  again  the  Coles'  function. 

The  most  straightforward  procedure  for  determining  the  heat 
transfer  coefficient  h  is  to  first  calculate  the  local  skin  friction 
coefficient  and  then  to  make  use  of  the  Reynolds  Analogy  in  a 
generalized  form.  The  latter  connects  the  Stanton  Number  with  the 
Reynolds  Number  by  means  of  the  expression 


S  « 


1  Cf 


P  c 


P 


“l  P  2 


(97) 


The  quantity  p  is  known  as  the  Reynolds  Analogy  factor  and  reflects 
the  difference  between  the  actual  value  of  S  and  that  determined  by 
Reynolds'  very  simple  analysis  in  which  he  demonstrated  that  p  =  1. 

In  general  p  is  primarily  a  function  of  Prandtl  Number  as  indicated 
by  the  example  of  Colburn's  analysis  of  heat  transfer  for  low  speed 
turbulent  boundary  layers  for  which  he  found  that 

P  =  (P)2/3  .  (98) 

This  of  course  gives  Reynolds'  result  when  P  =  1.  It  would  be 
reasonable  to  expect  that  in  the  case  of  a  high  speed  flow,  the 
Reynolds  Analogy  factor  would  depend  on  the  Mach  Number  and  on 
same  characteristic  temperature  ratio.  In  Harkness's  analysis  he 
showed  that 
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(99) 


p 

_ _ w_ 

p  ■  1  +  c  * 

where  c  ia  the  coefficient  defined  by  Eq.  (9l). 
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V.  EXPERIMENTAL  METHODS  FOR  DETERMINING  BOUNDARY  LAYER  CHARACTERISTICS 


A  variety  of  techniques  have  been  developed  for  the  experimental 
determination  of  boundary  layer  characteristics,  the  nature  of  which 
usually  depends  on  the  type  of  information  that  is  desired.  Probably 
the  earliest  of  these  methods  as  applied  to  flat  plates  was  repre¬ 
sented  by  the  work  of  Kernpf  [Ref.  (20)]  who  suspended  a  complete 
plate  from  the  carriage  of  a  towing  tank  and  measured  the  total  drag 
force  with  the  plate  parallel  to  the  direction  of  motion.  By  varying 
the  length  of  the  plate  and  the  towing  speed,  a  range  of  Reynolds 
Numbers  could  be  obtained.  When  applied  to  measurements  made  in  air, 
this  method  was  not  very  satisfactory  because  of  the  reduced  magnitude 
of  the  forces  to  be  determined. 

5.1  The  Momentum  Deficit  Method 

This  procedure  is  directly  related  to  the  momentum  integral 
theory  which  connects  the  frictional  drag  on  a  portion  of  the 
plate  surface  with  the  loss  of  momentum  within  the  boundary  layer 
up  to  the  point  where  the  measurements  are  made.  Since 


it  is  only  necessary  to  determine  the  local  velocity  u  as  a  function 

2 

of  y  and  to  compare  the  integrated  value  of  pu  dy  with  that  correspond¬ 
ing  to  the  free  stream  velocity.  In  the  case  of  a  compressible  fluid 
it  is  of  course  necessary  to  calculate  the  density  of  the  fluid  at 
each  point.  Values  of  the  velocity  are  readily  determined  by  means 


41 


of  a  Pitot  tube  attached  to  a  traversing  mechanism  so  as  to  give  u 
as  a  function  of  y.  In  the  incompressible  case  the  simple  Bernoulli 
Equation  provides  a  relation  between  Pitot  pressure  and  velocity, 
while  in  the  compressible  case,  variations  in  density  must  be  taken 
into  account,  including  the  use  of  the  Rayleigh-Pitot  formula  for 
supersonic  flow. 

In  order  to  minimize  errors  due  to  the  dimensions  of  the  Pitot 
tube,  the  visual  practice  is  to  maintain  the  transverse  dimension,  or 
diameter  if  the  tube  is  circular,  as  small  as  possible  consistent 
with  the  avoidance  of  excessive  time  lags  in  the  transmission  of 
changing  pressures  to  the  recording  instrument.  Such  considerations 
are  also  important  in  seeking  to  obtain  velocity  measurements  close 
to  the  plate  surface  and  if  possible  within  the  laminar  sublayer. 

For  these  reasons  seme  experimenters  have  utilized  Pitot  tubes  with 
a  flattened  cross  section,  the  minimum  dimension  obviously  being  in 
the  direction  normal  to  the  boundary  surface.  In  all  cases  care 
must  be  exercised  to  insure  that  the  Pitot  tube  is  properly  aligned 
with  the  main  flow  direction.  Since  the  measurements  are  being  made 
in  a  region  with  a  velocity  gradient,  there  may  be  seme  effect  on 
the  true  pressure  reading,  or  alternately  there  can  be  a  shift  in 
the  position  of  the  effective  c  '*nter  of  the  Pitot  tube  opening.  For 
tubes  of  reasonably  small  diameter,  however,  this  correction  is 
usually  negligibly  small. 

In  the  case  of  a  laminar  boundary  layer  or  within  the  laminar 
sublayer  of  a  turbulent  one,  it  may  be  desirable  to  obtain  readings 
at  a  sufficient  number  of  points  near  the  surface  so  as  to  detemine 
the  velocity  gradient  with  reasonable  accuracy.  The  shear  stress  at 
the  wall  being  given  by  the  relation 
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t  =  u 


(101) 


It  la  apparent  that  this  goal  may  be  difficult  to  attain,  since 
du/dy  must  be  determined  as  the  slope  of  an  experimentally  obtained 
curve  of  u  =  f(y) .  This  problem  is  avoided  in  the  case  of  pipe 
flow  by  relating  the  wall  shear  stress  to  the  longitudinal  pressure 
drop,  the  latter  being  a  quantity  which  can  be  measured  with 
reasonably  good  accuracy  by  means  of  static  pressure  taps  at  two 
stations  along  the  length  of  the  pipe. 

5.2  Local  Shear  Stress  Measurements 


In  view  of  the  fact  that  the  theoretical  calculations  of 
turbulent  boundary  layer  characteristics  directly  involve  the  wall 
shear  stress,  it  is  apparent  that  a  method  for  determining  the  value 
of  this  stress  has  significant  advantages.  The  integrated  values 
such  as  mean  skin  friction  coefficient  can  then  be  calculated  with 
good  accuracy.  Such  considerations  have  led  to  the  use  of  floating 
element  balances  in  which  a  portion  of  the  boundary  surface  is 
separated  from  the  surrounding  area  and  supported  by  means  of  a 
mechanism  through  which  the  shear  force  on  the  isolated  element  can 
be  measured.  Two  principal  types  of  shear  stress  balances  have  been 
developed,  the  initial  work  being  that  of  Liepmann  and  Dhawan 
[Ref.  (2l)]  who  used  a  null- reading  arrangement.  Since  some  clearance 
around  the  shear  element  is  required,  this  system  has  the  advantage 
that  the  element  is  always  in  the  same  position  when  the  reading  is 
taken.  The  other  type,  used  extensively  in  the  ARL  program 
[Ref.  (22)],  is  of  the  displacement  variety,  in  that  the  shear  element 
is  allowed  to  move  a  small  distance  in  the  flow  direction.  This 
element  is  supported  by  cantilever  springs  to  which  is  attached  the 
core  of  a  linear  variable  transformer.  Displacement  of  the  core 
causes  a  change  in  the  output  voltage  and  the  relation  between  this 
quantity  and  the  applied  shear  force  can  be  readily  determined  by 
calibration.  The  displacement  type  of  balance  is  particularly 
suitable  for  use  in  blow-down  wind  tunnels  where  short  running  times 
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are  involved.  With  both  types  of  balances,  it  is  essential  that  the 
shear  element  be  carefully  aligned  with  the  surrounding  surface.  A 
displacement  in  either  direction  normal  to  the  main  surface  can 
introduce  large  errors,  while  an  angular  misalignment  can  likewise 
be  extremely  troublesome .  Considerable  time  under  the  ARL  program 
has  been  devoted  to  a  study  of  these  alignment  problems,  with  the 
results  being  covered  primarily  in  a  report  by  O'Donnell  [Ref.  (23)] 
and  in  a  technical  journal  paper  by  Westkaemper  and  O'Donnell 
[Ref.  (24)].  While  not  directly  connected  with  the  Navy  supported 
program  at  ARL,  a  fairly  extensive  study  was  undertaken  for  the  NASA 
Langley  Research  Center  leading  to  the  design  and  construction  of 
improved  balances  to  be  operated  in  supersonic  flows  with  moderately 
high  rates  of  heat  transfer. 

In  general  the  desigu  of  a  satisfactory  skin  friction  balance 
requires  consideration  of  the  range  of  forces  to  be  measured,  these 
in  turn  being  dependent  primarily  on  the  Reynolds  Number  and  the 
Mach  Number.  Hence,  each  balance  must  be  treated  as  a  custom  built 
item  adapted  to  the  particular  flow  environment  with  which  the 
experiments  are  concerned.  It  is  always  desirable  to  have  as  high  a 
degree  of  sensitivity  in  the  balance  as  possible,  this  sensitivity 
being  controlled  by  the  size  of  the  shear  disk,  the  dimensions  of 
the  supporting  flexures,  and  the  response  of  the  transformer.  Since 
it  is  desired  to  obtain  shear  force  readings  effectively  at  the  center 
of  the  balance  disk,  a  reduction  in  its  size  is  always  desirable. 

This  reduction  is  limited,  however,  by  the  fact  that  the  forces  to 
be  measured  become  decreasingly  small,  while  at  the  same  time  the 
flexure  dimensions  become  so  small  as  to  introduce  severe  stability 
and  vibration  problems.  The  experience  at  ARL  has  shown  that  the 
minimum  disk  diameter  should  not  be  less  than  0.3  in.  and  preferably 
of  the  order  of  1.0  in. 
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5.3  The  Preston  Tube 

A  modification  of  the  Pitot  tube  application  was  developed  in 
the  early  years  of  aerodynamic  research  by  Sir  Thomas  Stanton  at  the 
British  National  Physical  Laboratory.  Stanton's  idea  was  to 
construct  a  total  pressure  probe  in  such  a  way  that  the  boundary 
surface  Itself  formed  the  inner  wall  of  the  tube.  The  overall 
dimensions  could  then  be  made  quite  small  so  that  the  tube  was 
capable  of  exploring  very  thin  boundary  layers  including  the  laminar 
sublayer.  An  extension  of  this  scheme  was  introduced  more  recently 
by  Dr.  J.  H.  Preston,  now  at  the  University  of  Liverpool  [Ref.  (25)]. 
In  his  application  an  ordinary  Pitot  tube  of  round  cross  section  is 
moved  toward  the  boundary  surface  until  it  is  actually  in  contact 
with  the  latter.  If  the  outside  diameter  of  the  nose  of  the  tube 
is  d,  then  the  center  is  at  a  distance  d/2  from  the  surface.  Assuming 
that  the  boundary  layer  is  turbulent  and  that  the  law  of  the  wall  is 
applicable,  Preston  then  reasoned  that  there  should  be  a  correlation 
between  the  surface  shear  stress  at  the  nose  of  the  tube  and  the 
pressure  reading  recorded  by  it.  As  a  part  of  the  work  reported  on 
in  Ref.  (ll),  Fenter  made  an  analysia  of  the  Preston  tube  in  a 
compressible  boundary  layer  flow,  extending  a  study  made  by  Hsu 
[Ref.  (26)]  for  the  incompressible  case,  which  includes  the  effects 
of  tube  diameter  and  wall  thickness.  The  significant  result  of  these 
analyses  is  that  the  tube  should  extend  well  into  the  region  in  which 
the  law  of  the  wall  is  valid  but  not  into  the  law  of  the  wake  region. 
Fenter' s  work  demonstrates  that  the  critical  or  maximum  diameter  of 
the  tube  is  given  approximately  as 


d 

cr 


~  0.226x  r  "y5 
~  t+1  x 


(102) 


Here  x  is  the  distance  from  the  plate  leading  edge  to  the  measuring 
station,  R  is  the  corresponding  Reynolds  Number,  and  t  is  the  wall 
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thickness  of  the  nose  of  the  tube.  As  was  the  case  with  the 
floating  element  balance,  the  AKL  Investigations  of  the  performance 
of  the  Preston  tube  also  included  the  effects  of  pressure  gradients 
and  heat  transfer  [Refs.  (27)  and  (28)]. 


VI.  RESULTS  OF  SKIN  FRICTION  AND  HEAT  TRANSFER  MEASUREMENTS 


In  general  the  experimental  phase  of  the  ARL  program  followed 
In  parallel  with  the  analytical  work,  with  the  former  covering  first 
the  use  of  the  momentum  deficit  method  on  smooth  flat  plates  under 
adiabatic  conditions.  These  investigations  were  next  extended  to 
measurements  on  plates  with  uniformly  distributed  roughness,  with 
both  the  floating  element  balance  and  the  Prestou  tube  being 
introduced  into  the  program  at  this  time.  The  consideration  of  more 
regular  types  of  roughness,  particularly  that  consisting  of  transverse 
Veegrooves,  was  studied  next  with  the  effects  of  heat  transfer  also 
Included,  with  the  smooth,  adiabatic  plate  being  included  as  a  limiting 
case  of  zero  heat  transfer  and  zero  roughness.  The  results  of  each 
of  these  various  phases  of  the  experimental  program  are  summarized 
in  the  sections  which  follow. 

6.1  Skin  Friction  Measurements  on  the  Smooth  Flat  Plate 


The  measurements  on  smooth,  adiabatic  plates  were  actually  an 
integral  part  of  the  program  reported  on  by  Wilson  in  Ref.  (7)  and 
involved  the  application  of  the  momentum  deficit  method  using  a 
traversing  Pitot  tube  across  the  boundary  layer  at  several  stations 
along  the  plate  center  line.  The  results  thus  led  to  determinations 
of  mean  skin  friction  coefficients  as  functions  of  Reynolds  Number. 
Tests  were  made  at  several  different  Mach  Numbers  in  the  low  super¬ 
sonic  range.  The  range  of  Reynolds  Numbers  extended  from  2  X  10^ 
to  19  x  10^,  with  corrections  being  made  for  the  position  of  the 
effective  leading  edge  of  the  plate  as  influenced  by  the  presence  of 
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a  turbulence  tripper  strip  near  the  physical  leading  edge.  The 
range  of  Mach  Numbers  extended  from  =  1.579  to  2.471,  these 
values  being  detemined  by  the  nozzle  blocks  available  for  the 
Ordnance  Aerophysics  Laboratory  wind  tunnel  in  which  the  tests 
were  run.  A  set  of  blocks  for  =  2.75  was  also  available  but  was 
not  used  in  the  boundary  layer  program  because  of  the  relatively  poor 
quality  of  the  test  section  flow. 

The  results  of  these  tests  can  be  presented  graphically  in 
condensed  foxm  by  making  vise  of  the  analysis  which  led  to  Eq.  (52). 
Thus  Pig.  14  is  a  plot  of  X  as  a  function  of  Y  with  the  latter  two 
quantities  defined  by  Fqs.  (49)  and  (50),  respectively.  According 
to  Wilson's  analysis,  this  plot  should  result  in  a  single  curve  which 
is  independent  of  Mach  Number.  The  experimental  values  for  each  of 
the  test  values  of  Mach  Number  are  indicated  separately  in  Fig.  14 
and  vexy  satisfactorily  confirm  the  analysis  which  resulted  in 
Eq.  (52). 

For  engineering  purposes  and  to  minimize  the  amount  of 
numerical  computation  required,  the  set  of  charts  developed  by 
Clutter  and  presented  in  Ref.  (9)  are  probably  to  be  preferred, 
since  they  yield  directly  the  values  of  Cp  for  specified  values  of 
the  Reynolds  Number,  with  a  separate  chart  being  available  for  each 
Mach  Number  between  M^  =  0  (incompressible  flow)  and  =  5.  These 
charts  also  include  the  effects  of  sand- grain  type  roughness, 
represented  by  the  ratio  of  the  plate  length  f,  to  the  mean  sand 
particle  dimension  k. 

Wilson's  measurements  of  velocity  across  the  boundary  layer 
were  studied  to  only  a  limited  extent  in  comparison  with  the  velocity 
profiles  predicted  by  the  analysis.  There  were  indications  that 
pointed  up  the  need  for  improved  assumptions  regarding  the  shear 
stress  and  mixing  length  distributions,  but  these  points  will  not  be 
discussed  in  detail  in  the  present  report  since  they  do  not 
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significantly  alter  the  relationship  between  skin  friction 
coefficient  and  Reynolds  and  Mach  Numbers. 

6.2  Skin  Friction  Measurements  on  Roughened  Plates 

The  initial  phase  of  the  ARL  investigation  of  roughness  effects 
begem  as  a  natural  extension  of  the  smooth  plate  experiments,  following 
lines  comparable  to  those  employed  by  Nikuradse  in  his  studies  of 
turbulent  flow  in  pipes.  The  sand- grain  type  of  surface  was  formed  on 
the  smooth  plate  by  first  coating  it  with  a  mixture  of  clear  varnish 
and  drier.  The  roughness  itself  was  formed  by  applying  grinding 
compound  of  various  sizes  to  the  plate  with  a  flocking  gun  with  care 
being  taken  to  obtain  as  uniform  a  surface  as  possible.  The  results 
of  these  initial  experiments  on  roughness  are  presented  in  a  report 
by  Shutts  and  Fenter  [Ref.  (29)],  the  tests  being  made  in  the  QAL  wind 
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tunnel.  The  Reynolds  Number  range  was  from  4  x  10  to  2  x  10,  while 
the  Mach  Number  varied  from  1.62  to  2.50  in  increments  of  approximately 
0.25.  Several  different  sizes  of  grinding  compound  were  utilized, 
the  pertinent  data  being  summarized  in  Table  III. 

A  typical  plot  of  mean  skin  friction  coefficient  as  a  function 

of  Reynolds  Number  is  shown  for  one  Mach  Number,  M^  =  2.00,  in 

Fig.  15.  Three  values  of  the  roughness  parameter  are  included  in 

this  plot  corresponding  to  values  of  R  -  0,  the  smooth  plate,  and 

x  x  r 

for  =  2.59  x  lw  and  4.01  X  10  .  The  theoretical  curves  shown 
in  Fig.  15  were  based  on  Wilson's  analysis  of  Ref.  (7)  in  the  case 
of  the  smooth  surface,  and  on  Fenter' s  treatment  of  uniformly  rough 
surfaces  as  presented  in  Ref.  (ll). 

The  test  program  represented  by  Ref.  (29)  included  boundary 
layer  velocity  surveys,  as  well  as  skin  friction  balance  determinations 
of  local  shear  stress.  The  latter  values  were  integrated  with  respect 
to  x,  the  distance  along  the  plate,  to  obtair  the  values  of  the  mean 
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skin  friction  coefficients  shewn  in  Fig.  15.  A  limited  comparison 
of  local  skin  friction  coefficients  with  theory  was  presented,  with 
agreement  comparable  to  that  obtained  for  the  mean  values. 

Since  these  tests  represented  the  first  experience  at  ARL  with 
roughened  surfaces  and  with  the  use  of  the  skin  friction  balance, 
it  is  considered  significant  to  note  some  of  the  special  problems 
resulting  from  the  introductions  of  the  grit-type  of  roughness  and 
the  balance  measurements.  The  grit  particles,  when  observed  under 
a  microscope,  were  extremely  irregular  in  shape,  so  that  a  statistical 
approach  was  indicated  to  obtain  the  mean  roughness  dimension.  For 
a  given  grit  as  many  as  250  samples  were  observed  microscopically 
and  the  average  dimension  of  these  samples  used  as  a  means  of 
determining  the  value  of  r,  the  mean  roughness  dimension.  In  the 
case  of  the  skin  friction  balance  measurements  on  the  roughened 
surface,  it  was  necessary  to  use  extra  care  to  insure  that  roughness 
particles  did  not  enter  into  the  gap  between  the  balance  disk  and 
the  surrounding  plate  area.  Such  fouling  would  of  course  prevent 
the  disk  from  experiencing  the  noimal  displacement  to  be  expected 
from  the  action  of  the  shear  force. 

Additional  phases  of  the  roughness  program  at  ARL  included 
extensions  of  the  work  reported  on  in  Ref.  (29),  with  particular 
attention  being  given  to  a  comparison  of  the  skin  friction  balance 
and  Preston  tube  methods  of  deteiminlng  local  values  of  the  skin 
friction  coefficient.  A  special  investigation  was  initiated  to 
determine  the  equivalent  roughness  of  surfaces  having  other  than 
the  grit  type  coating,  which  was  accomplished  by  Installing  suitably 
machined  Inserts  in  a  basic  flat  plate  model.  These  tests  were 
conducted  in  both  the  QAL  and  ARL  wind  tunnels.  Another  aspect  of 
the  program  involved  measurements  on  a  cone- cylinder  model  in  the 
supersonic  wind  tunnel  of  the  Vought  Aeronautics  Division  of 
Ling-Temco-Vought,  Inc.  The  details  of  these  various  segments  of 
the  program  were  presented  in  a  series  of  ARL  reports,  but  the 
essential  findings  Eire  well  summarized  in  Ref.  (ll). 
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6.3  Skin  Friction  and  Heat  Transfer  Measurements  oa  Roughened  Plates 


The  final  phases  of  the  AHL  roughness  effects  program  were 
extended  to  include  measurements  under  heat  transfer  conditions.  At 
the  same  time  it  was  decided  to  standardize  on  a  Veegroove  type  of 
roughness  as  being  representative  of  machined  surfaces  resulting 
from  actual  production  procedures .  The  initial  work  along  these 
lines  is  covered  in  the  report  by  Young  [Ref.  (12)]  already  cited 
in  connection  with  the  analytical  studies.  His  experiments  involved 
the  addition  to  the  plate  of  a  coating  of  tin- lead  solder,  followed 
by  a  rolling  process  in  which  Veegrooves  of  the  desired  dimensions 
were  impressed  upon  this  coating.  Since  the  basic  plate  model  was 
fabricated  of  copper,  considerable  difficulty  was  experienced  in 
attempting  to  machine  such  grooves  into  the  plate  itself.  The 
rolling  process  applied  to  the  solder  coating  appeared  to  avoid 
the  galling  that  occurred  with  the  machining  procedure  as  Indicated 
by  microscopic  examination  of  the  completed  surface. 

The  particular  pattern  of  Veegrooves  selected  for  this  study 
had  a  90  deg  angle  at  the  peak,  so  that  the  height  of  the  projection 
was  always  one- half  of  the  base  width.  The  longitudinal  axis  of 
the  grooves  was  perpendicular  to  the  plate  center  line  and  the  flow 
direction.  While  seme  degree  of  temperature  elevation  could  be 
introduced  into  the  wind  tunnel  air  supply,  provision  was  also 
made  for  internal  cooling  of  the  plate  model  so  as  to  increase  the 
range  of  the  ratio  of  wall  to  free  stream  temperature.  The  internal 
cooling  system  was  also  designed  so  as  to  control  the  plate  to  a 
uniform  temperature  over  its  entire  surface.  The  specific 
dimensions  used  in  this  series  of  experiments  were  roughness  heights 
of  0.005,  0.010,  and  0.030  in. 
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The  instruaentation  employed  in  the  program  had  to  meet  the 
requirement  of  simultaneous  measurements  of  local  skin  friction  and 
of  heat  transfer  rate.  This  objective  was  accomplished  by  providing 
two  holes  in  the  plate  at  a  distance  of  12.5  in.  aft  of  the  leading 
edge  with  centers  located  1.0  in.  on  each  side  of  the  longitudinal 
center  line.  One  of  the  two  holes  was  designed  to  accommodate  a 
floating  element  balance,  while  the  other  provided  for  the 
installation  of  a  plug- type  calorimeter.  Both  the  balance  shear 
disk  and  the  calorimeter  plug  had  diameters  of  1.00  in.  The  nominal 
Mach  Number  of  the  wind  tunnel  air  flow  was  4.90,  while  the  tempera¬ 
ture  ratio,  T^,  ranged  from  about  5.2  to  2.9.  Variations  in  this 
ratio  could  be  obtained  by  adjusting  the  stagnation  temperature  of 
the  tunnel  air  flow,  the  wall  or  surface  temperature  of  the  plate, 
or  both.  Heat  transfer  rates  were  determined  by  a  transient  method 
in  which  the  surface  temperature  of  the  calorimeter  plug  was  decreased 
abruptly  by  injection  of  about  10  cm^  of  cold  water  on  its  surface. 

The  injection  tube  was  then  pulled  up  out  of  the  main  flow  and  the 
calorimeter  and  surrounding  plate  temperatures  were  recorded  at 
frequent  time  Intervals .  The  value  of  the  Stanton  Number  was 
calculated  on  the  basis  of  the  temperature- time  gradient  at  the 
Instant  when  the  calorimeter  and  plate  temperatures  were  the  same. 

In  this  manner  any  errors  due  to  incomplete  insulation  of  the 
calorimeter  disk  were  minimized.  The  details  of  this  procedure  are 
fully  described  in  Ref.  (12). 

The  skin  friction  balance  readings  were  taken  simultaneously 
with  the  calorimeter  readings,  but  the  Pitot  surveys  across  the 
boundary  layer  had  to  be  made  in  separate  runs.  It  should  also  be 
noted  that  the  temperature  ratios  for  the  runs  with  the  different 
degrees  of  roughness  were  not  exactly  the  same  throughout  the 
series,  due  to  difficulties  in  maintaining  precise  temperature 
control.  Since  all  other  parameters  such  as  Reynolds  Number,  Mach 
Number,  and  degree  of  roughness  were  essentially  constant,  it  would 
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have  been  preferable  to  Interpret  the  data  on  the  basis  of 
variations  in  cf  and  c^  with  temperature  ratio.  Actually  the  test 
program  was  executed  with  a  constant  value  of  the  plate  temperature 
at  =  555  °R  and  the  temperature  ratio  varied  by  adjusting  the 
tunnel  stagnation  temperature  and  the  ambient  air  temperature,  T^. 

This  procedure  led  to  variations  in  Reynolds  Number  with  temperature 
ratio  for  a  given  measurement  station,  although  this  variation  was 
not  large. 

A  consideration  of  all  of  these  factors  finally  brought  the 
investigators  to  the  conclusion  that  it  was  best  to  consider  the 
local  skin  friction  coefficient  and  the  Stanton  Number  as  functions 
of  the  roughness  Reynolds  Number,  represented  by 

Rr  “  u*  k^/2  .  (105) 

In  addition  to  the  data  given  in  Ref.  (12),  involving  Veegroove 
heights  up  to  0.050  in.,  results  obtained  by  Mann  [Ref.  ( 50) ]  are 
also  Included  in  these  plots,  which  are  shown  here  in  Figs.  l6  and  17. 
The  initial  work  of  Young  was  found  to  carry  the  roughness  height 
approximately  to  the  so-called  threshold  value,  and  Mann's  study 
was  therefore  undertaken  in  order  to  extend  the  roughness  heights  to 
larger  values,  specifically  0.060  in.  and  0.090  in. 

The  plot  of  Fig.  l6  shows  the  variation  in  local  skin  friction 

coefficient,  cf,  with  the  roughness  Reynolds  Number  Rr.  Although  not 

strictly  correct,  it  is  assumed  that  the  temperature  ratio  is  constant 

for  each  of  the  curves  drawn.  If  the  values  of  T  /T,  are  considered 

w'  1 

as  "standard"  for  the  smooth  plate  runs,  the  percentage  variation  in 
temperature  ratio  in  general  does  not  exceed  ±5$,  although  there  are 
one  or  two  cases  of  roughened  plate  where  the  variation  was  as  high 
as  7  to  12$.  Another  factor  involved  in  interpreting  these  results 
is  based  on  the  fact  that  the  so-called  "smooth  plate"  was 
undoubtedly  not  completely  smooth  in  an  aerodynamic  sense.  Thus  while 
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the  test  values  of  are  plotted  as  corresponding  to  R  =0,  the 
horizontal  lines  drawn  through  these  points  serve  to  indicate  that 
they  should  be  displaced  slightly  to  the  left.  Profilograph 
measurements  over  a  representative  area  of  the  smooth  plate  would 
have  served  to  establish  an  appropriate  value  of  kr,  but  such 
equipment  was  not  available  for  this  purpose  during  the  test 
program.  A  similar  comment  also  applies  to  the  heat  transfer  data 
shown  in  Fig.  17  in  which  the  Stanton  Number  S  is  plotted  against  Rr» 

In  addition  to  the  small  variations  in  temperature  ratio  between 
one  surface  and  another,  a  more  significant  factor  should  be  noted, 
i.e.,  that  there  is  a  substantial  variation  in  Reynolds  Number  between 
runs.  Thus  for  the  adiabatic  wall  condition,  the  nominal  Reynolds 
Number  is  R  =  14  x  10  ,  while  for  the  runs  at  a  temperature  ratio 
T^/P  =2.8,  the  Reynolds  Number  has  decreased  to  5  X  10^.  Thus, 

Figs.  l6  and  17  are  really  two  parameter  families  of  curves  involving 
both  temperature  ratio  and  Reynolds  Number  as  parameters.  A  more 
useful  correlation  of  these  measurements  would  have  been  obtained  if 
the  Reynolds  Number  had  been  maintained  at  a  constant  value  while 
the  temperature  ratio  was  varied  for  a  given  degree  of  roughness. 

This  would  presumably  require  that  the  tunnel  temperature  be  held 
constant  and  the  temperature  ratio  be  varied  by  means  of  the  model 
cooling  system.  Such  an  approach  might  have  resulted  in  severe 
restrictions  on  the  range  of  temperature  ratios  that  could  be 
covered.  The  execution  of  such  a  test  program  in  a  continuous- flow 
rather  than  a  blow-down  wind  tunnel  would  also  be  preferable,  since 
once  a  stable  flow  was  established  in  the  tunnel,  the  temperature 
ratio  could  then  be  varied  systematically  during  a  single  run 
without  interruption. 

Using  the  data  shown  in  Figs.  l6  and  17,  values  of  the 
Reynolds  Analogy  factor,  p,  as  defined  by  Eq.  (97)#  were  computed 
and  plotted  as  functions  of  temperature  ratio  in  Fig.  18.  While 
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rough  mean  curves  could  be  drawn  through  the  data  points  for  each  of 
the  different  roughness  values,  the  extreme  scatter  of  the  data  made 
It  undesirable  to  do  this.  Also  such  curves  seem  to  exhibit  no 
logical  or  systematic  trends  with  Increasing  roughness.  The  remedy 
for  this  situation  Is  undoubtedly  to  be  found  In  the  more  elaborate 
test  procedure  outlined  above,  particularly  one  Involving  more  pre¬ 
cise  control  of  all  flew  parameters  as  well  as  the  plate  temperature. 
Similar  results  are  given  in  a  technical  Journal  paper  by  Young  and 
Westkaemper  [Ref.  (3l)]  but  based  only  on  the  Initial  group  of 
roughness  values  reported  on  in  Ref.  (12). 

As  mentioned  earlier  In  connection  with  the  outline  of  Young's 
research  program,  his  experimental  work  also  Included  Pitot  probe 
surveys  of  the  velocity  distribution  across  the  boundary  layer. 

These  results  were  expressed  initially  in  terms  of  the  velocity 
ratio,  u/u^,  as  a  function  of  the  normal  distance  y  from  the  plate 
surface.  Curves  for  the  smooth  plate  and  the  different  degrees  of 
roughness  are  shown  in  Pig.  19  for  the  case  of  the  Insulated  plate, 
that  is,  Tw  =  or  T =  5.2.  The  data  for  the  smooth  plate 
and  for  kf  /  0  are  so  close  together  that  only  the  curve  for  the 
smooth  case  is  shown.  The  other  cases  show  a  steepening  of  the 
velocity  gradient  near  the  wall  with  Increasing  roughness  which 
of  course  Is  consistent  with  the  expected  increase  In  skin  friction 
shear  stress.  A  second  plot  for  T^/T^  =  3.8  is  shown  In  Fig.  20 
and  Indicates  little  change  due  to  the  change  in  thermal  conditions. 
In  both  cases  the  free  stream  Mach  Number  was  the  usual  value  of 
=  4.93.  These  velocity  measurements  may  be  readily  transformed 
into  the  usual  semilogarithmic  plots  with  u/u*  plotted  as  a  function 
of  u*y/2;  such  a  plot  for  T^/T.^  =  3.8  is  shown  in  Figs.  20  and  21 
for  the  two  cases  just  mentioned.  In  this  form  the  effect  of 
roughness  is  more  readily  apparent  but  the  data  are  so  irregular 
in  character  that  no  firm  conclusions  can  be  drawn  from  them  as  to 
the  precise  magnitude  of  the  change.  Again  it  would  appear  that 
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there  Is  a  real  need  for  experiments  In  a  continuous- flew  wind  tunnel 
under  more  precise  control  of  temperature.  One  would  hope  that  such 
experiments,  when  repeated,  would  lead  to  quantitative  data  on  rough¬ 
ness  effects  comparable  to  Rlkuradse's  work  on  low  speed  pipe  flow. 

Although  this  report  was  prepared  primarily  as  a  summary  of 
the  ARL  contributions  to  the  turbulent  boundary  problem,  including 
roughness  and  heat  transfer  effects,  an  effort  has  been  made  through¬ 
out  to  Introduce  references  to  other  significant  work  In  this  area 
of  recent  date.  In  this  connection  the  present  report  might  be 
closed  with  a  reference  to  a  Symposium  on  Compressible  Turbulent 
Boundary  Layers  held  at  the  NASA  Langley  Research  Center  In 
December  1968  [Ref.  (32)].  The  Proceedings  of  this  meeting  provide 
a  state-of-the-art  stannary  for  ccmpressible  flow  comparable  In  many 
respects  to  the  earlier  Stanford  report  on  incompressible  flow.  In 
both  cases  there  appears  to  be  a  vigorous  battle  being  waged 
between  two  schools  of  thought,  first,  those  who  favor  the  continued 
use  of  Integral  methods  for  the  prediction  of  skin  friction  and  heat 
transfer  rates,  and  second,  those  who  advocate  a  return  to  a  more 
fundamental  approach  which  might  eventually  lead  to  more  rational 
theories.  The  large  amount  of  empiricism  involved  In  the  first 
method  is  one  of  the  major  factors  in  limiting  the  usefulness  of 
the  integral  approach.  Carefully  planned  and  conducted  experiments 
on  turbulent  boundary  layers  In  all  aspects  of  the  speed  spectrum 
would  do  much  to  provide  a  sounder  physical  basis  for  analysis  than 
that  now  available. 
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TABLE  I 

ROUGHNESS  FUNCTION  PAR i 
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5 
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— 
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TABLE  II 


THE  NUMERICAL  REPRESENTATION  OP  THE  LAW  OF  THE  WALL  AND 
THE  LAW  OP  THE  WAKE  AS  DEVELOPED  BY  COLES  [REF.  (3)] 
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FIGURE  4 

LAW  OF  THE  WALL  FOR  COMPRESSIBLE  FLOW 


FIGURE  5 

MACH  NUMBER  PARAMETERS  FOR  COMPRESSIBLE  FLOW 
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WILSON’S  THEORY  -  r  =  1.00 
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FIGURE  12 

VARIATION  IN  ROUGHNESS  FUNCTION 
WITH  ROUGHNESS  REYNOLDS  NUMBER 
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